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ABSTRACT 


The oceanographic community is moving towards unmanned 
autonomous vehicles to gather data and monitor scientific 
sites. The mission duration of these vehicles is dependent 
primarily on the power consumption of the propulsion system, 
the control system and the sensor packages. 

A customized propulsion thruster is designed. This 
includes a specialized propeller tailored to ABE and a matched 
motor and transmission. A non-linear lumped parameter model 
of the thruster is developed and experimentally verified. The 
model is used to predict thruster performance and compare the 
design thruster with other variants of propeller and 
motor/transm1ission combinations. 


The results showed that there is a trade-off between 
rapid dynamic response and power conservation. For the 
typical ABE trajectory, the designed thruster provides good 
dynamic response and the lowest power consumption of all the 
modelled thruster units. 
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Chapter 1 


Introduction 


1.1 Motivation 


The Woods Hole Oceanographic Institution (WHOI) is 
currently developing an autonomous underwater vehicle (AUV) 
for scientific survey of the ocean floor. This vehicle has 
been designated ABE, for the Autonomous Benthic Explorer. 

As the oceanographic community explores the ocean floor, 
the reliance on manned submarines to maintain ongoing 
experiments has become restrictive. ALVIN, WHOI's manned 
deep submersible, 1s unable to undertake new research, due to 
an exhaustive schedule maintaining experiments that are 
already in progress. ABE 1s being developed in order to free 
up assets such as ALVIN and JASON, WHOI’s unmanned 
submersible, by performing routine data collection and 
surveying at remote ocean bottom scientific sites. 
Additionally, since ABE can operate without a nearby surface 
ship operating as a tender/control ship, ABE will also free up 
valuable research vessel time. In these respects, ABE will 
complement the existing capabilities of tethered and manned submersibles. 


JR 


ABE will be a long endurance vehicle with a typical 
deployment of one year in length. This compares to an on- 
station endurance of eight hours for a manned research 
submersible (ALVIN) and around a month for tethered robots 
(JASON). During the deployment, ABE will observe a relatively 
small area ( on the order of square kilometers ) of the ocean 
floor at frequent intervals. Figures 1.1 and 1.2 show the 


preliminary ABE configuration and a typical mission profile. 


During a mission, ABE will remain in a semi-dormant state 
for the majority of the time. At a predetermined interval or 
in response to a trigger event, ABE will wake, and preform a 
photographic survey along a preprogrammed flight path. Upon 
completion of the survey, ABE will return to its mooring, 


power down and wait for the next survey time. 


Black & White Vertical Thrusters 
Cameras 


Flotation 
Pods 


~ 
` z ‘S 
` f 4 
4 p 4 
elena 


Docking 
Mechanism 


Side Thrusters 33 
Strobes 


Color Camera Instrument & 
Battery Housing 





Figure 1.1 The Basic ABE Design. [1] 
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In order for ABE to have sufficient battery power for a 
one year mission, careful attention must be paid to developing 
a highly efficient propulsion system. This research 1s 


motivated by the need to develop this propulsion system. 


Autonomous 
Benthic Explorer 


* self powered 

* pre-programmed 

* long life (1 year) 

* can move 50 km at 1 knot 

e carries films 4 CCD camera 

* camies CTD & transmissometer 


Glide unti over target ; 
© 3 * survey may be event triggered 


—— — = 
<p 


+ (3) Spiral down over 


<P bottom beacon 
When 20m above © + 
bottom, move away 


Ti Perform regular stereo photo survey 
of a pre-programmed region 
Or survey region around | 
another area ot Interest 


HP 


Latch onto beacon “hitching post (6) XE 
and switch to low power mode 





Figure 1.2 The Typical ABE mission profile 
will consist of four phases: Descent (1-6), 
Sleep (6), Survey (7-8), and Ascent. [1] 


1.2 Research Objectives 


There are two major components of a AUV propulsion 
system. The first is the mechanical/hydrodynamic system 
commonly referred to as a thruster or propulsor. The thruster 
consists of an electric motor, a transmission and a 
propeller/duct. The second component is the electronics and 


the algorithm used to control the mechanical system. 


to 


The thruster must be optimized for the specific vehicle. 
This involves designing a specialized propeller and duct 
assembly then matching a suitable motor and gearbox to take 
advantage of the motor's high speed efficiency and the 
propeller's low speed efficiency. This is done to maximize 
the conversion of electrical energy to thrust. 

The control system must then be designed to provide the 
best possible dynamic response, while taking advantage of the 
mechanical system's most efficient operating conditions. 

The objective of this research is to develop a propulsion 
system that provides good dynamic response while maintaining 


the high efficiency needed by ABE to perform 1ts mission. 


1.3 Outline of Thesis 


Chapter 2 presents the design of an efficient propeller 
for the ABE vehicle. Chapter 3 examines the selection of a 
motor and gearbox matched to the ABE propeller. In Chapter 4, 
a lumped parameter model of the thruster is developed. 
Chapter 5 contains the experimental verification of the model. 

Chapter 6 compares the steady state performance of the 
designed ABE thruster with other thruster units. Chapter 7 
examines the trade-offs that must be made between power 
efficiency and beneficial dynamics. Chapter 8 summarizes the 
results of the thesis and provides recommendations for further 


research. 
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Chapter 2 


The Propeller and Duct 


The goal of this chapter is to develop a propeller and 
duct combination for ABE that optimize the power conversion of 
the thruster. 

2.1 Introduction to Propeller Design 

The propeller design presented here is based on the 
Massachusetts Institute of Technology's (MIT's) propeller 
development software. This software uses Lifting Line Theory 
and optimum circulation to solve the complex hydrodynamics of 
a marine propeller. [2] 

The process used to design the propeller for ABE is as 
follows: 

1. Calculate the drag of the tentative ABE 


vehicle at the desired operating velocity. 


2. Determine the thrust required. 


3. Determine the physical constraints of the 


vehicle that effect propeller size. 


4. Estimate the wake field near the propeller. 
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5. Perform a parametric study using MIT's 
Propeller Lifting Line program (MIT-PLL). 


Minimize Power and torque for a given thrust. 


6. Use MIT’s Propeller Blade Design program 
(MIT-PBD10) to calculate the blade shape from 


MIT-PLL's optimum output. 


2.2 ABE Vehicle Constraints 


2.2.1 ABE Physical Considerations 


The ABE vehicle consists of two buoyancy pods supporting 


a single instrument cylinder underneath. The buoyancy pods 


are twenty one inch diameter series 58 bodies, and the 


instrument cylinder 1s a twelve inch diameter streamline 


cylinder. All three bodies are seven feet long and they are 
connected by a series of struts. There are seven 
thruster/propulsor units on ABE. Three are main propulsors 


mounted at the stern of each cylinder. The remaining four are 
thrusters for attitude and depth control. Two of these last 
four thrusters are mounted vertically, and the other two are 
mounted athwartships. ABE is designed to have a cruising 
speed of one knot, and a minimum sprint capability of two 
knots. The power budget allows 100 watts for normal 


propulsion with a peak of 200 watts available. The propulsors 
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are limited to a diameter of eighteen inches by the space 


avallable on ABE. 


2.2.2 Vehicle Drag 


The buoyancy pods are made up of twenty-one inch, serles 
58, streamlined shapes. Each pod will have a surface area of 
PS LEE? The dragon this Eyperor body tor relatively slow 
speeds is primarily due to friction on the surface. The drag 


can be calculated from the following equation: 


Drag = Ž pU? CS (2,1) 


Where the drag coefficient, C,, 1S approximately equal to the 
flat plate frictional drag coefficient, C,;. For ABE at one 
knot in seawater, the Reynold’s Number is =6 X 10° and the 
corresponding Cs is $0.007.[3] S is the wetted surface area, 
U is the velocity through the water, and p is the density of 
the seawater. Using this equation, the drag of each of the 
buoyancy pods at one knot is 0.56 lbf. The cylindrical 
instrument case will be fitted with streamlined nose and tail 
cones and equation (2.1) holds for this case as well. The 
drag for the twelve inch diameter instrument case at one knot 
Isf 005447 IDE. The combined drag of the three main body 


sections Of ABE is 1-0- IDF: 


L7 


Attention must be paid to the struts connecting the body 
sections. Careful design of these structures Can reduce 
vehicle drag significantly. If it is assumed that the struts 
are well designed and have a streamlined cross-section, then 
equation 2.1 above can be used to calculate the drag. The 
term S is now equal to the entire wetted surface of the strut 
(top and bottom surfaces). Assuming nine feet of struts with 
one inch maximum thickness and five inches width, the 
calculated drag is =0.15 lbf. 


However, 1f the struts are poorly designed the drag 


increases by an order of magnitude. Assuming a l-inch 
circular cross-section, the drag is now calculated by: 
"lt 
Dragsrror = > PU? dC, (2.2) 


Where the drag coefficient, C,=1.0, for laminar flow and d is 
the diameter of the strut. For nine feet of circular struts 
the drag 1S approximately 2.25 lbf. 

For ABE, the struts are assumed to be made up of 
streamlined shapes for two-thirds of the total length and 
bluff sections for the remaining length. The bluff sections 
account for joints, fouling, and imperfections in the strut 
sections. Using this distribution the drag due to the struts 
is 0.8 Ibf armik knot: 

A large portion of the drag will come from the attitude 


control thrusters. The. “attitude Controll thrusters are 
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comprised of four thrusters mounted perpendicular to the flow 
streamlines when ABE is traveling in the forward direction. 
These thrusters will be used to make minor corrections to the 
vehicle's depth and heading. These thrusters will be ducted 
to prevent fouling and impact damage. Each unit will have a 
3 inch diameter, streamlined motor case and a 3 inch by 18 
inch duct that will be streamlined for forward motion. Using 
the above equations, each thruster will have a drag of 0.2 
St 

The total vehicle drag for a forward speed of 1 knot is 
3.2 1bf. For the purposes of propeller design, each thruster 
will be designed to give 3 lbf at 1 knot. This allows ABE to 
maintain a full speed capability in the event of the failure 
of one of the three main propulsion thrusters. Additionally, 
ABE will have an excess propulsive force available during 


normal cruising conditions. 


2.3 Optimum Efficiency 


Using the thrust of 3 lbf determined from the drag, the 
vehicle’s speed of one knot, an estimated shaft speed of 100 
rpm, and the 1.5 foot propeller diameter, the ideal efficiency 
can be determined from the Kramer diagram, Figure 2.1 [2]. 
The entering arguments for the diagram are A (a form of the 


ship advance coefficient, J;) and C,, the thrust coefficient. 
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For ABE these coefficients are: 
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Figure 2.1 Kramer Ideal Propeller Efficiency 
[2] 
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The ideal efficiency for a three blade propeller 
operating under these conditions is n;=83%.[2] This 
represents the maximum efficiency possible, neglecting viscous 


forces. 


2.4 A Parametric Study using MIT-PLL 


Using the information already discussed, we are almost 
ready to begin using MIT-PLL to start developing some chord 
and thickness distributions for the ABE propeller blade. In 
order to enter MIT-PLL, initial chord and thickness 
distributions must be assumed and a inflow velocity field 
needs to be determined. 

The initial chord and thickness distribution were assumed 
to be linear. At the hub, the chord 1s 2 inches, the 
thickness is % inch. At the blade tip, the chord is ¥% inch 
with a tenth of an inch thickness. 

For an accurate inflow field to be determined, extensive 
model testing must be performed. Since this type of testing 
1s beyond the scope of this paper, a simple inflow field is 
assumed. This field is shown in Figure 2.2. 

With this information, a parametric study of several 
potential propellers was conducted. The parameters used in 
this study are: (1) Number of Blades, (2) RPM, and (3) Blade 
Shape (in a qualitative sense). Thes goal SWECO fine eine 


maximum efficiency and the minimum power required to generate 
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Wake Inflow Profile 


0.6 0.65 0.7 0.75 





Fraction of Free-Stream Velocity 


Figure 2.2 The assumed inflow velocity 
field 


the required 3 lbf thrust, with an acceptable blade shape. 


The first step is to determine the number of blades for 


the propeller. To do this, the required thrust was varied 
from the 3 lbf value calculated above. This determined two 
things: 1) How dependent efficiency calculations are on 


thrust, and 2) how dependent these calculations are on the 
number of blades. Table 2.1 shows the results of these runs. 
For each run, the optimum RPM was determined, the effective 
wake recalculated, circulation optimized and chord lengths 
optimized. The output was run through the cycle until 
horsepower, torque and efficiency remained constant. Figure 
2.3 . shows the plot of efficiency versus thrust for V,=1 knot 
and various numbers of blades. The figure shows that the 
number of blades has little effect on the efficiency over 


ABE's expected operating range. 
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Thrust_ li 0.5 
Bld: 3 








| 


Torque 





Table 5 No Tunnel Used ın Calculations) 
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Over the range of thrust from 1 to 6 1btf, the three, foun 
and five blade propellers show only a slight variation from 
each other. This difference can be considered statistically 
insignificant. The three, four and five blade propellers were 
used in the next stage of development. All propellers have 


comparable results. 


Efficiency for Various Blade Numbers 


fee oad ee 


5 Blades 


3 Blades 


3 
Thrust, lbf 





Figure 2.3 Efficiency, N versus required Thrust. 


The evolution of the final blade took many iterations. 
The first 3-blade propeller had a predicted efficiency of 
=70%. (This was calculated without a tunnel around the 
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propeller.) The resulting chord distribution had an extremely 
narrow blade with a very sharp tip. This blade was not 
physically suited for use in the marine environment. Any 
fouling or impact with an obstruction (eg. fish) would have 
destroyed the blade. 

Several chord distributions were used. By taking the 
output of one run, the planer blade shape was plotted (based 
on Ne: chord distribution). From this plot, a new, more 
rugged chord distribution was developed. Thirsa Chord 
distribution was used for the next MIT-PLL run as the initial 
bil'zde input. For cach run of MIT-PLL, the thrust (3 1bÉ),. the 
diameter (1.5 ft) and the inflow velocity field were held 
constant. The RPM, circulation, and chord lengths were 
optimized. After several iterations, a final 3-blade 
propeller with a 66.1% efficiency anda 4-blade propeller with 
a 66.4% efficiency were chosen. These propellers showed the 
highest efficiency, with durable blade dimensions. 

The three blade propeller is presented since a 
commercially available propeller similar to the designed 
propeller was readily obtained. The commercial propeller is 


presented in section 2.7 


2.5 MIT-PBD10, Blade Shaping 


After using MIT-PLL to determine the desired circulation, 


chord and thickness distributions, the actual blade shape and 
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camber that will develop the desired circulation must be 
calculated. The code used for this is MIT’s Propeller Blade 
Design program, PBD10. To enter MIT-PBD10, a rake and skew of 
the blade is required in addition to the MIT-PLL output. 
Since rake has negligible effects and serves no purpose for 
the ABE vehicle, no rake is used. The primary purpose of skew 
is to balance (by phase shift) unsteady forces on the 
propeller. Since the forces on ABE are small and the speeds 
of operation are low, the unsteady forces are neglected anda 
small amount of skew (8° at the tip) was added to aid in 
obstruction shedding. 

The recommended default values for MIT-PBD10 were used 
initially. These values determine the nature and extent of 
the wake field. The two dimensional blade cross-section shape 
was chosen as a NACA a=0.8 mean line. Slight modifications 
were made to the extent and contraction of the wake field in 
order for PBD10 to run smoothly in this particular case. The 
PBD10 output for K, K, and the induced velocities are similar 
to what was described by PLL. Figures 2.4 through 2.7 show 
the resulting blade shape determined by the above process. 
Figures 2.8 and 2.9 show a planer view of the blade, with and 


without skew. 
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The summary of the propeller designed 1s: 


Number of Blades: Z=3 


Hub Radius: Rog 
Tip Radius: R,=8 .4" 
Duct OD: D=18" 
ID: d=17" 
Optimum RPM: N=94.5 


al knot & 3 1bf 


Horsepower: P=.008 
Torque: t=0.46 ft-lbf 
Efficiency: n=0.661 
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Figure 2.4 Chord Distribution 





Camber (Fo/C) Distribution 
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Figure 2.5 Camber Distribution 


Blade Thickness Distribution 
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Figure 2.6 Thickness Distribution 
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Blade Pitch Distnbuuon 


Pitch, Inches 


Planer Blade Shape, No Skew 


“A 
Q 
ie 
E 
nv 
2 
G 
3 
x 
V 
2 
= 
CQ 


-2 0 2 


Blade Width, Inches 





Figure 2.8 Planer Blade Shape, No Skew. 


Planer Blade Shape, with Skew 
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Blade Width, Inches 
Figure 2.9 Planer Blade Shape, Light Skew. 


2.6 The Duct 


The next step in the ABE propulsor development is to 
design the duct. The primary purpose for the ducts on ABE is 
to prevent fouling of the propellers and the shafts during 
prolonged deployments. The goal is to design a low drag duct 
that has a minimal effect on the vehicle’s hydrodynamics and 
propulsion efficiency. The MIT-PLL program traces the neutral 
streamlines at the tip of the propeller. If a duct is placed 
along these streamlines, the duct will have no effect on the 
propeller. If a slight angle of attack is placed on the duct, 
lift can be achieved on the duct section. This lift, when 
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summed around the entire duct, becomes a propulsive force (as 
opposed to drag). Using MIT T-PLL; seewovec tml he war ssligat 
propulsive force was designed. The neutral nose-tail angle of 
attack is 4.4 degrees. The -angle “Or attack. Lor S547 duct 
propulsion is 3.8 degrees. The cross-sectional shape for the 
duct is the NACA 0008 Basic thickness form. [4] Figures 2.10 


and 2.11 show this duct. 


E 


5% Duct Shape 





Figure 2.10 5% Propulsive Duct 


Duct Cross-section 
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Figure 2.11 Duct Cross-Sectional Shape 
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2.7 The Experimental Propeller 

The test propeller used in collecting data was a Michigan 
Wheel 18-inch diameter, 16-inch pitch, Sailer ™, 3-Blade 
propeller. This propeller was chosen because it most closely 
matches the chord and pitch distribution of the designed 


propeller. Figures 2.12 through 2.14 describe this propeller. 


Expenmental Chord Distnbution 
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Figure 2.12 Experimental Propeller Chord 
Distr but ron 


DÒ 


Experimental Blade 
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Figure 2.13 Experimental Propeller Pitch 
Dictribution 





Experimental Blade 
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Figure 2.14 Experimental Propeller Planer Blade Shape 
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Chapter 3 


The Motor and Gearbox 


3.1 ABE Vehicle Constraints 

There are several serious constraints placed on the 
mechanical design of the motor and transmission by ABE's 
mission environment. The two primary considerations are the 
ambient pressure of the ocean bottom and the limited supply of 
electrical power provided by onboard batteries. 

In order to survive the extreme pressure (10,000 psi), 
the motor assembly must elther be pressure tolerant or 
isolated from the pressure. In order to isolate the motor, a 
heavy pressure vessel must be constructed and the motor sealed 
inside. There must be a shaft seal around the output shaft. 
This seal must be leak proof at an extremely large 
differential pressure. Such a seal 1S expensive and produces 
a large added load on the motor. This load reduces system 
efficiency ina dramatic way. Pressure tolerant motors also 
require shaft seals. However, since the internal pressure of 


the motor casing is maintained at a few pounds above ambient 
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pressure, the differential pressure (DP) is very low. For low 
DPs, the added load is very low. Additionally, seals that 
operate with high differential pressure are prone to slight 
leakage. Any leakage of seawater over an extended period of 
time can result in failure of electric devices. For these 
reasons, shaft seals with high differential pressure should be 
avoided in this application, and a pressure tolerant motor 
should be used. 

The disadvantage of the pressure tolerant motor is an 
increased loss of power due to windage. The windage comes 
from the fluidic drag on the rotor of the motor due to the 
presence of a fill fluid. The £111 fluid is a non-conductive 
fluid maintained at a few psi above ambient pressure and it 
surrounds and fills the motor. This prevents the highly 
conductive and corrosive seawater from entering the motor. 

Since ABE is a battery powered vehicle with a bus voltage 
of 48 VDC, a DC motor is the obvious choice for the prime 
mover of the thruster. 

The transmission must be selected to match the motor, 
which is most efficient at high speeds, to the propeller, 
which has high efficiency at lower speeds. The gearbox must 
also be of sufficiently high quality and precision to minimize 


the losses due to the gearing. 
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3.2 DC Brushless Motors 


3.2.1 Brushed versus Brushless Motors 

DC motors come in two main types: Brushed and Brushless. 
Brushed motors are the most common DC motor. They have a 
mechanical arrangement of split rings and brushes called a 
commutator. The commutator switches the voltage applied to 
the coils depending on the position of the rotor. This 
switching keeps a positive force on the rotor to keep it 
rotating. A brushless motor relies on electronics to provide 
the proper commutation to the motor based on feedback from an 
external rotor position detector. 

In a high pressure environment, spring loaded brushes 
experience increased wear and a tendency to hydroplane on the 
non-conductive fill fluid. The hydro-planing leads to brush 
chatter and an increased heat load due to the increased 
electrical resistance. Brushed motors have a short life 
expectancy in the high pressure environments. 

Since the commutation on brushless motors is accomplished 
electronically, they do not suffer from any of the above 
problems. They are, however, much more expensive and require 
complicated (and expensive) controllers. Due to the 
importance of longevity and reliability in ABE's mission 


environment, DC brushless motors will be used. 
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3.2.2 Motor Equations 

In order to evaluate the motor, the efficiency of the 
motor is calculated. The motor chosen for ABE is a Pittman 
elcom © 5100 series DC brushless motor. This motor was chosen 
for size and rated capacity. A schematic of the motor, 
defining the variables used in the equations, is shown in 


Pagure 3... 


Controller V, 





Figure 3.1 Motor Schematic 


The motor constants and parameters for this motor are: 


Torque Constant: 


A 
+3 
Il 


0.173 Nm/Amp 


Back EMEF Constant: K: 


0,173 V/(rad/sec) 


Coil Resistance: RO 4.85 ohms [oa 


DIS 


For steady state operation, the following equations describe 


the operation of a DC motor: 


Back EMP: Me 
Motor Torque: US 
Motor Current: OR A ne 


Where V, is the supply voltage and @, 1s the angular velocity 








of the rotor. Efficiency can be calculated as: 
n, = Pout = Tn®m (3.1) 
Lia Val y 


From equation 3.1, it can be seen that efficiency 
increases with motor speed. This neglects the effect of 


windage which increase with speed. 


3.3 The Fill Fluid and Windage 

The fluid used for compensation 1s Halocarbon © 0.8 cSt 
mead. "Tharsis a silicon ol based fluid with a viscosity 202 
less than water. Of the commercially available, non- 
conductive fluids, this fluid provides pressure compensation 
and has the lowest viscosity. For the purposes of efficiency 
Galicilation, “the fluid flow around ‘the rotor will be 
considered to be laminar flow between two parallel surfaces. 
Brom testing of ¡other motors. with this. fluid, (a -linear 
relationship exists between torque due to drag and angular 
velocity. The proportionality constant in this relationship 
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is K, = 5.175 x 107 Nm/(rad/sec). The torque lost (t,) and the 


power lost (P,) due to windage are: [6] 


and P, = K o 


t, = Kio 3 > 


wW wo m 


3.4 The Gearbox 

Gearboxes are available in a wide variety of types, gear- 
ratios and efficiency. Gearbox efficiency is highly dependent 
on the manufacturer’s tolerances and construction procedures. 
For the Pittman motors, a variety of planetary gearboxes are 
available. These gearboxes will be used to determine the 
desired gear ratio for the ABE thruster. In order ue 
calculate the desired gear ratio, descriptive equations of the 
performance of the gearbox must be determined. The equations 
must be written in terms of gear ratio. 

To formulate the equations it is assumed that a planetary 
gearbox 1S made up of an arbitrarily small stage. A certain 
number of these stages are stacked in order to get the desired 
reduction ratio. Each stage has a specific gear reduction 
(1.1:1) and a specific efficiency (n;). The complete gearbox 
then has a gear reduction ratio of (1.1)":1 and an efficiency 
of n; where n is the number of stages (including fractional 
stages) needed to get the desired reduction. For the Pittman 
gearboxes, two advertised gearboxes are a 4:1 and a 17.33:1 
gear ratio with efficiencies of 80% and 64% respectively. If 
the 4:1 gearbox is used as a baseline, n=14.545 and n;=0.985. 
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If these numbers are used to calculate the efficiency for the 
aa: Lo géearbox, n-=23:93 and tie etrtiecieney is "033, since 
this 1s a good match, the following equations will be used to 


describe the gearbox: 


CO A O la (3.3) 


3.5 Matching the Motor with the Propeller 

In order to evaluate a proper match between the propeller 
and the motor, a set of equations describing the motor, 
windage and gearbox must be evaluated for various gear ratios. 
Combining the optimum propeller RPM and torque from chapter 2 
with the component equations developed above, a system can be 


evaluated. Figure 3.2 shows the block diagram of the system. 


Propeller 


Tau, (LOSS) 
w 





Figure 3.2 Thruster System Block Diagram 
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The descriptive equations are: 





Torop = 0.46 £tIDÎ = 0.62386 Nm 
= = rad 
W crop = 94.5 RPM = 9.896 oe 
DE prap reno 
i> NeNc 1 7” *"G™ prop 
Um "a Ti £ Ty Ty = KyO n 
T 
O, = di Iy n K. 
Vs = Vaacx + RrIm Vaack = Kz@y 


N; = (1.1)" 6 = (.985)2 


Volt 
rad 
( sec ) 


Nm 
K, = 0.173—— K 2210 7 3 
T Amp E 





Nm 


(2) 


SO LO 


These equations reduce to: 


Va AO 
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T 
O A O 
”  (1.1)”Kp(0.985)” A 
with the component efficiencies defined by the following: 


Motor efficiency without windage: 


n. = Power out _ bem 
m Power in Vai 


Motor efficiency with windage: 


n = Power out = Cm dm us Ty 
De: Power in We 


Gearbox efficiency: 


n; = (0.985)? 


Total system efficiency: 


Nr = Nw G 


These efficiencies were calculated for various gear ratios. 
Figure 3.3 shows a plot of these efficiencies and Figure 3.4 
Shows an enlargement of the total efficiency curve. By 


inspection of this figure, the ideal gear ratio is about 10:1. 
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Efficiency 


Total Eîfisieney ion... Motor Effic 


Gear Ratio 


15 


Gear Reduction 


Figure 3.4 Total System Efficiency 
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3.6 Sensitivity of the Analysis 

This analysis was carried out using several different 
stage reduction ratios and efficiencies. The effects of 
changing these parameters are summarized below. 

1. Changing the stage reduction ratio only from 1.1:1 to 
4:1 had little effect on the gear ratio where the peak 
efficiency occurred. It did effect the height of the peak 


Significantly. 


2. Varying the stage efficiency alone (fixed 1.1:1 stage 
reduction ratio) had little effect on the peak location while 
the stage efficiency was above 95%. Below 95% the peak moved 


to lower gear ratios. 


3. Varying both the stage reduction ratio and efficiency 
in a coordinated manner to maintain 80% efficiency at a 4:1 
gear ratio, had no significant effect on both peak location 


and height. 


In general, the model of the gearbox and motor developed 
in this chapter is relatively insensitive to the assumptions 
made about the gearbox. The windage loss dominates the peak 
efficiency curve, provided a gearbox with sufficiently high 


efficiency is used. 
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Chapter 4 
The Lumped Parameter Model 


of The Propulsor 


4.1 Introduction 

ABE will incorporate a complex control system designed to 
allow the vehicle to operate independently for periods of up 
to twelve months. In order to accomplish this mission, ABE 
must have extremely efficient thrusters and an efficient 
control algorithm. The control algorithm must be based on a 
simple, yet accurate model of the thruster. In this chapter 
a lumped parameter model of the thruster 1s developed and 
presented in bond graph notation. A simulation of this model 
was performed using MATLAB, and the results are compared to 


experimental data in Chapter 5. 


4.2 Description of the Thruster Unit 


The thruster under consideration, consists of a DC 
brushless motor, a controller for the motor, a 10:1 reducing 
gearbox, a three blade propeller mounted ina duct. Figure 
4.1 shows a schematic diagram of this assembly. 
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Propeller 
& Duct 


Tau ,,, (LOSS) Tau, (LOSS) 





Figure 4.1. Schematic representation of the 
thruster unit. 


4.2.1 The Propeller and Duct Assembly 


The propeller and duct assembly as a unit, has the single 
most significant impact on the overall efficiency of the 
thruster. Therefore, the model of this unit determines the 
accuracy and utility of the overall model. Figure 4.2 shows 
the schematic representation of the propeller and duct. In 
order to describe the hydrodynamic relationships, the 


following simplifying assumptions are made [7]: 


de. The ambient fluid is inviscid, incompressible and of 


constant density. 


oe The gravity effects on the fluid are negligible. 
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ce The flow at the inlet and outlet is parallel to the 


Lhrister axis. (Irrota Cioni 

de The only energy storage in the fluid is by kinetic 
energy. 

5a The kinetic energy of the ambient fluid is 
negligible. 

6. The ambient pressure is P, and it acts equally at 


the inlet and the outlet. 





Figure 4.2. The propeller and duct 
assembly, Schematic Representation. 


In order to model the hydrodynamics, first, consider the 
Kinetic Co-energy, Ex . This 1s the typical "Physics Book 


Kinetic Energy" { % mv? }. 
Be= (ev)( 2) (4.1) 


Where pV is the mass of the fluid enclosed by the duct. 


However, this volume must be corrected to take account of the 
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added mass effect of the surrounding fluid. This correction 
involves a concept known as Added Mass or Virtual Mass. Added 
mass 1S a phenomenon that occurs when a body (or a fluid 
element) moves through a fluid. An additional inertia effect 
is added to account for the effort required to "push" the 
fluid out of the way of the moving body. For the case of a 
body moving slowly through a static fluid, the added mass 1s 
equal to the volume of the body times the density of the fluid 
(ie. the mass of the displaced fluid) [8] and [9]. For this 
model (where the density is constant) this is done by setting 
V to twice the actual enclosed volume. Q is the volumetric 
flow rate. A is the cross-sectional area of the duct. Now, 
define the pressure momentum, I, as: 
la la). E (4.2) 

“Note that there is a linear relationship between T and 
O. This is analogous to the standard translational definition 
of momentum, p=mv, and leads to the numerical equality of the 


Kinetic Energy, E, and the Kinetic Co-energy: 


Ez = frao = ¿PY E) (4.3) 








£ far i, ar? _— a [pvoĵ _ pvo _ x 
E = fod = f su __ 2pV 2pV| A? Ane È 
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In Summary: 





=p: - PON. 221? al 
ENEE ee T (4.5) 
r= £9 And Qi al (4.6) 
A? PV 


Now consider a power balance for the propeller and duct. 


Power in: 


The power flow into the duct comes from three major 
components: 1) The driving motor, 2) Any Kinetic Energy 
flowing into the thruster inlet, and 3) The velocity/ 
opposing-force product at the inlet ( work done by the fluid 


entering). 


1. The power input from the motor/gearbox is the product 


of the torque, T, and the angular velocity, 0. 


2. Since it is assumed that the ambient fluid is at roS 


there is no kinetic energy flow into the thruster. 


3. The velocity/force power is the product of ambient 


pressure P, and the cross-sectional area and the fluid 


velocity at the inlet. 
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n 


PAT (P.a)( 5 eG OSO (4.7) 


Power out: 


The power flows out of the duct by two processes. The 
first is the force/velocity work at the discharge. The second 
is the kinetic energy of the discharged fluid. This latter 


flow is called the convected kinetic energy. 


1. Since the cross-sectional area is constant and the 
average velocity of the fluid is the same at the inlet and 
outlet, this term is the same as the corresponding term at the 


Inlet a PO: 


2. The convected kinetic energy is the kinetic energy per 


Unit volume times. Che volumetric -Elow rate: 


E 272 
K = Convected Kinetic Energy = o| = ha (4.8) 
p 


P - pọ + 2T lol (4.9) 


out O > p y? 


Note: the absolute value preserves the sign of the convected 


kinetic energy to allow for flow reversal. 


Dal 


Balance of Power: 


The net rate of change of the energy in the fluid 


contained 1m the duct as equal te Pe. 


one Ol IO) 
LE, = Pin- P = on Peroni Sd 


did 2 p V? 2p V? 
Since the only method of energy storage within the duct 


1s through kinetic energy: 





Ean E a) =A Tiro - AIO (4.11) 
pV Poe 


ps - Hel (4.12) 


However, 0=0(F)=0(@%) therefore we can write FT as I(@) and 


reduce equation (12) as follows: 


1. Define p as the pitch of the propeller. This quantity 
is also known as the advance of the propeller. Specifically, 
pis the distance the propeller travels axially per revolution 
in an ideal est lw: 

2. Define N, the propeller's efficiency, is equal to 1-0; 
where 6 1S the slip of the propeller: [7] 


az 


ODA = 250 


o â (ALS) 
w pA 
Now we can write [I and O in terms of 0. 
g= 2pna (4.14) 
2 T 
where Q = (rev per sec)*(pitch)*(efficiency)* (Area) 
= Volumetric Flow rate 
r- 2negvo o) 
2TA 
Using equations (4.6) and (4.15), we can rewrite equation 
(412) as: 








p- 27r _ ATI) _ 2nt _ polo) (4.16) 


pA 2p V? pA 2A 


This equation can easily be written in terms of O: 


2 
dra An T _pnaolo| i7) 


By py * amy 


From these equations we get the bond graph shown in Figure 


4.3. 


SS 





Figure 4.3. Bond Graph for Propeller-Duct Assembly. 


The Equations for the Propeller and duct are: 














2 Tec 
[ee Sn - Pon (4.18) 
R is Defined such that Pon > ae (4.19) 
AT V 
E 7 ee 4.20 
°° pv Le eli | 
w =- 222 T=- 2% (4.21) 
pA ppA 


The dissipator (R) is often referred to as a "Bernoulli 
Resistor". In this model, the convected kinetic energy is 
dissipated to the ambient fluid. By assuming that the ambient 
fluid remains at rest, it is assumed that the ambient fluid 
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acts as an infinite sink for this energy, much the same as 1s 
considered for the thermal energy dissipated by an electrical 
e iStor anlo.: 

The parameters for the propeller are: p=0.41 meters/rev., 
Propeller efficiency is 60%, A= 0.164 m?, corrected volume is 


02033 mn 


4.2.2 The Gearbox 


The motor is connected to the propeller shaft by means of 
a gearbox. This gearbox has a planetary gear arrangement 
yielding a 10:1 reduction of speed. This gearbox is necessary 
to increase overall efficiency since the motor operates most 
efficiently at high speeds and the propeller operates best at 


lower speeds. The gearbox arrangement 1s shown in Figure 4.4. 





= 


Tau, 


Omega Omega 


Figure 4.4. Gearbox Schematic Diagram 
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The gearbox used has an 80% efficiency at top speed 
(input shaft maximum speed is 10,000 RPM). This means that 
there is a dissipative element as well as a torque 
multiplication. The inertial load of the gearbox will be 
included with the motor’s rotor inertia. Defining T, as the 
motor output torque, T, as the gearbox output torque, @ as the 
input shaft speed, and @ as the output shaft speed, the 
gearbox can be modelled as shown in the block diagram of 


Figure 4.5 
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Figure 4.5. Bona Graph for Gearbox. 


The equations describing the gearbox are: 





0.27 
R¿= m N,=10 (4.22) 
W max 
0.27 W 
CEN: tn T vg = ——— = 0.2Tt, (4.23) 
W max 


S6 


Jg will be combined with other inertias into an equivalent 


inertia. 


4.2.3 The Motor 


The motor used is a Pittman elcom © DC Brushless motor. 
It 1s controlled using a +/- 10 VDC control voltage and a 
Supply "voltage ranging from 30 to 80 Volts Dc. In “Chars 
application the supply voltage is chosen to be 48 VDC. In 
order to simplify the model, the supply and control voltages 
are "tied" together. This idealizes the controller portion of 
the model by assuming no losses in the amplification. The 


motor diagram 1s shown in Figure 4.6 


48 Vdc Supply 


DC Brushless 
Controller 


Motor 


Control Voltage O to »/- 10 Vadac 





Figure 4.6. The Motor Assembly. 
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In order for the motor to operate on the ocean floor, e 
must be surrounded by a non-conductive fluid that is 


compensated to be at the same pressure as the environment 


(10,000 psi). The fluid in use is a silicon based oil called 
Halocarbon © 0.8 cSt. This oil adds a significant windage 
loss to the motor ( = S 2x 10 na Nmn aa aea 


The motor is then modelled by Gyrator, a mechanical 
inertia and dissipator, and an electrical resistance and 
inductance.{11] The gyrator constant is 0.173 Nm/Amp or 07 we 
V/(rad/sec). The stator resistance R,=4.85 ©, the stator 
inductance L,= 1.65 mH, the combined rotor and gearbox moment 
of inertia J,=38.6 x 10 7° Kg m2. 

These numbers are from the manufacturer's specifications. The 


resulting bond graph 1s shown in Figure 4.7 
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Figure 4.7. The motor bond graph. 
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The resulting equations model the motor: 


Um n Ti n J Ò n 7 Kym 
Ti A KrIm 
Vback 3 KeWa 
E 
Iy = TL Veontrot - Vback ~ RrIy] 


4.3 The Complete Thruster Model 


The complete model is shown in Figure 4.8 
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Figure 4.8 
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(4.24) 


(4.25) 


(4.26) 


(4.27) 


In summary the equations are: 
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The output of the system is thrust. Thrust is equal to 
the time rate of change of momentum. For our system the 
thrust is equal to the convected momentum, specifically, the 


product of momentum per unit volume and the volumetric flow 


racos 
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3 | | 
Thrust = ( momentum per unit volume ES, = 70 as = Di 
p 


Thrust can also be written in terms of oO: 


Apn polv 
An? 


TDrusSt = 


The bond graph representation of Figure 4.8 can be put 


into block diagram form. The resulting block diagram of the 


complete thruster unit is shown in Figure 4.9. 






Desired 
Thrust 


>) 


Figure 4.9 The Complete Model Block Diagram 
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This model contains the electrical dynamics of the motor 
as well as the hydrodynamics of the propeller. The electric 
time constant L.,/R, =0.34 milliseconds. This is extremely fast 
compared to the hydrodynamic time constant which is on the 
order of a second. Therefore the hydrodynamics dominate the 
electrical and the later can be neglected without degrading 
the model. The resulting model is shown in Figure 4.10 . The 
model presented here has been simplified to correspond to the 


model used in reference [12]. 


Desired 
Thrust 





Figure 4.10 The Simplified Model Block Diagram 


Where Rh, Js meg lected cinee kh < e anda 
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c, = AP°N°p a= PUPA 
An 167 Toy 


and Torque = B vV, and Power = Tw = BV.W where B is a motor 


constant with units of Newton Meters per Volt. 


4.4 Simulation 

The model of Figure 4.10 was used in MathWorks© MATLAB™ 
Program. This program was used to simulate the response of 
ee “thruster to” various “anputs. The results of these 


Simulations are compared with experimental data in Chapter 5. 


4.5 Summary 

In this chapter, a lumped parameter model was developed 
to describe a underwater vehicle thruster unit. From this 
lumped parameter model, a simplified block diagram model was 
derived. This model will be used to predict the response of 
an actual thruster unit under test at the Woods Hole 
Oceanographic Institution. The three parameters of the final 
model ( C,, A, and B) will be calculated and then tuned to get 
the best match between the actual response and the simulated 


response. 
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Chapter 5 


Experimental Verification 


The next consideration is to show that the model 
developed in Chapter 4 is competent to describe the actual 
thruster. This will be done in three steps. First, by 
equating the steady state responses of the model and the 
experimental thruster. Second, the model will be tuned to 
match the thrusters step response. Finally, the actual and 
predicted responses to several frequencies of sinusoid will be 
examined. The model’s ability to correctly predict the 


sinusoidal response will be evaluated. 


5.1 The Experimental Setup 


For the purposes of model verification, an MFM ™ DC 
brushless motor (electrically comparable to the Pittman motor 
described in Chapter 4) with a 10:1 gearbox, was mounted ina 
housing filled with Halocarbon fluid. The experimental 
propeller (Michigan Wheel 3-blade 16-inch pitch propeller) was 
mounted on the gearbox output shaft and supported radially by 
an external journal bearing. A frame was constructed to hold 


the motor housing and bearing assembly. The frame is 
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supported at a single pivot point with a lever arm extending 
upward to allow force measurement to occur out of the water. 
Figure 5.1 shows the arrangement of the frame and the force 
measurement sensor. The force sensor is an S-type load cell 


rated at 500 pounds-force. 


Load Cell 





Figure 5.1 The Experimental Test Stand 


After initial testing, the test frame was reinforced to 
reduce the oscillatory effect of the vertical lever arm 
compliance. This reinforcement significantly reduced the 
"ringing" of the thrust measurement. In order to remove 
gravity effects, the test assembly was balanced before each 
data run, so the vertical lever arm was straight up and down. 
The force sensor was arranged to measure the horizontal 
component of force at the end of the vertical arm. This 
arrangement decouples all gravitational effects of the motor 


and propeller from the force measurement. 
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5.2 Steady State Response 


The constants of the model developed in Chapter 4, (C,, 
a, and B) were adjusted so the model accurately predicts the 
correct steady state response of the thruster for thrust, 
angular velocity, and power input required. This process 


determined C, and PB directly. 


5.3 Step Response 


The step response data collected with the experimental 
thruster was used to determine the QA parameter. The modelled 
response to three separate step inputs is shown in Figure 5.2 


and the actual response is shown in Figure 5.3 (a,b,&c). 


Desired and Actual Thrust 


l ES 





Time, Seconds 


Figure 5.2 Model Step Response 
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(a) 76 Newton Step 





(b) 55 Newton Step 





(c) 5 Newton Step 


10 Squares 


0.5 Second 
3 Actual Thruster Step Response 


1 Tick Mark 
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The model and the actual thruster were given sinusoidal 


inputs of five different frequencies between 0.1 Hz and 1.6 


The model response is shown in Figures 5.4 through 5.8 


HZ. 
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sub-plot (a). The actual response is shown in subplot (c) of 
the same figures. Note that the actual response is inverted. 


The time scale is .1 second per finest square. The force scale 


is 5.7 N per smallest division. 
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SS Sp a SS ES SS 





Time Scale: 10 Squares = 1.0 Seconds or 1 Tick Mark = 0.5 Seconds 





Figure 5.4 Sinusoidal Input T=10 seconds. (a) Model Input and 
Response, (b) Actual Input, (c) Actual Response 
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Time Scale: 10 Squares = 1.0 Seconds or 1 Tick Mark = 0.5 Seconds 





Figure 5.5 Sinusoidal Response, T=5.0 Seconds, (a) Model Input 
ana Response, (b).Actual Input, (ec) Actual Response 
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1.0 Seconds or 1 Tick Mark = 0.5 Seconds 


Figure 5.6 Sinusoidal Response, T= 2.5 Seconds, (a) Model 
Input and Response, (b) Actual Input, 


(c) Actual Response 
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Time Scale: 10 Squares = 1.0 Seconds or 1 Tick Mark = 0.5 Seconds 





Figure 5.7 Sinusoidal Response, T= 1.25, (a) Model Input and 
Response, (b) Actual Input, (c) Actual Response 
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Figure 5.8 Sinusoidal Response T=0.625, (a) Model Input and 
Response, (b) Actual Input, (c) Actual Response. 
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By inspection of the sinusoidal responses, it is seen 
that the model accurately predicts the magnitude and the time 
lag of the actual thruster. This lag is firmly rooted in the 
hydrodynamics of the propeller. The time scale of this lag 1s 
the time necessary to develop the helical wake field down 
stream of the propeller after changes occur that effect the 
propeller's angular velocity. The model predicts the 
magnitude of this lag within 0.1 seconds of the actual 
response. The magnitude prediction is degraded for 
frequencies over 1 Hz where it under predicts the actual 


thrust developed by up to 40%. 


5.4 Summary 


The model developed in Chapter 4, and tuned in this 
chapter to match steady state and step responses, does a good 
job at predicting the dynamic response of the thruster. 
During the experimental studies, the actual thruster response 
is highly dependent on the adjustment of the controller to 
account for the inertia of the propeller. The controller used 
has two modes of operation: Torque control and Velocity 
control. The velocity control mode was very sensitive to the 
propeller inertia. If the controller is not well tuned, the 
velocity control mode becomes unstable. This risks damage to 
the mechanical linkages and uses a large amount of electrical 


power. Since this undesirable trait was not readily 
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correctable, all the comparisons use the controller's torque 
control mode. The model was tuned to best fit the controller 
when the controller was matched to the experimental propeller. 
In the next chapter, the model and controller were tuned to 
other motor and propeller combinations. The method of this 


chapter was used to tune all configurations tested. 
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Chapter 6 
Steady State Performance Comparison 


of Several Thruster Units 


6.1 Introduction 


The purpose of this thesis is to design the best thruster 
unit for ABE. In order to verify the performance, the 
experimental thruster unit will be compared to other 
combinations of propellers and motor/gearbox units. In this 
chapter, four separate propellers and two motor/gearbox units 
are tested to determine steady state Power required to obtain 
a certain thrust at a bollard pull. The test setup described 
in Chapter 5 was used to gather the data. The results of 
these tests are compared to determine which thruster has the 
lowest power consumption over the desired range of thrust. 

The motors/gearbox units under test were: 

1. MFM Technology Inc. Series SM64 DC Brushless 
M@bor witirta 102) «Gearbox: 
2. Pittman elcom Series 5100 DC Brushless 


Motor with a 4:1 Gearbox. 
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The Propellers used in these tests were: 
1. JASON Thruster Propeller; a 2-blade, 10- 


inch diameter propeller with a 4-inch pitch. 


2. 18-inch diameter, 2-blade, model airplane 


propeller with a 6-inch pitch. 


3. 18-inch diameter, 2-blade, model airplane 


propeller with an 8-inch pitch. 


4. 18-inch diameter, 3-blade, Michigan Wheel 
Sailer Marine propeller with a 16-inch pitch. 
(This 1s the EXPERIMENTAL Propeller) 
Each combination of these components was tested in the 
two controller modes, Torque and Velocity control. 
The “experimental" thruster consists of the MFM motor 


with the 10:1 gearbox and the Michigan Wheel propeller. 


6.2 The ABE Thruster and Other Thruster Units 
The thruster units will be categorized by the propeller 
used. The discussion starts with the smallest propeller and 


ends with the experimental propeller. 
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6.2.1 The JASON Propeller 


The first propeller considered is the propeller from a 
JASON vehicle thruster. This propeller is used by the 
tethered underwater vehicle JASON and is equivalent to a 
standard 2-blade Mercury outboard motor propeller. Figure 6.1 
shows the response of the JASON propeller to each 
motor/gearbox unit and controller mode. In this chapter the 
graph labels are decoded as follows: 

The first letter indicates the motor unit under 

test, M corresponds to the MFM motor (10:1 

reduction) and P corresponds to the Pittman Motor 


(ae yr eduction): 


The middle 2 or 3 letters/numbers indicate the 
propeller used. JAS = the Jason Propeller, 86 = the 
18" diameter 6" pitch airplane propeller, 88 = the 
S*ibiuch* dxrplane propeller, 316 = the Michigan 


Wheel 3-blade 16" pitch propeller. 


The last letter indicates the controller mode, T = 


torque control, R = velocity control. 


Examination of Figure 6.1 shows for the range of -10 to 
+10 lbf thrust, the Pittman motor/ 4:1 gearbox in torque mode 


outperforms the other methods of powering the JASON propeller. 


1 


For the 100 Watt ABE power budget limit, this unit provides a 


range of thrust from -8 lbf to +9 lbf at bollard. 


Propeller JASON 





Figure 6.1 The Jason Propeller Response 


6.2.2 The 86 Airplane Propeller 


The second propeller is the 18-inch diameter model 
airplane propeller with a 6-inch pitch. This propeller is 
extremely asymmetric and is designed to be operated in only 
one direction. Figure 6.2 shows the steady state data for 
this propeller. Note the large difference in power required 
between forward and reverse thrust. For the +/- 10 lbf thrust 
range all thruster unit combinations are comparable, with the 
MFM motor in torque mode being the best. For ABE’s 100 Watt 


limit, this propeller can provide ™—S or ee 12, IDE 
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Figure 6.2 The 86 Airplane Propeller 
Response 





6.2.3 The 88 Airplane Propeller 


The 18-inch diameter model airplane propeller with an 8- 
inch pitch 1s very similar to the 86 propeller above. Figure 
6.3 shows the data collected for this propeller. Again the 
MFM motor in torque mode outperforms the other combinations. 
The range of thrust available within the 100 Watt limit is -7 
lbf to 13 lbf. This is slightly better than the 86 propeller 


performance. 
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Figure 6.3 The 88 Airplane Propeller 
Response 





6.2.4 The Experimental Propeller 


The last propeller considered is the experimental 
propeller described in Chapter 2. This propeller, when 
matched with the designed motor/gearbox unit, outperforms the 
4:1 gearbox in both controller modes. The range of thrust 
avallable under the 100 Watt limit is -7.5 lbf to 11 1bf. 
This propeller has the best bi-directional response of any of 
the tested propellers. Figure 6.4 shows the data for this 


propeller. 


6.3 Steady State Comparison 


Due to the limitation of the controller mentioned at the 


end of Chapter 5, the velocity control mode was not used for 
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Figure 6.4 The Experimental Propeller 
Response 
the final comparisons of the thruster units. The velocity 


mode shows a marked instability that is highly dependent on 
the inertia of the attached propeller. Thise instability 
causes an excess of power to be drawn at various speeds of 
operation, including stopped. Since the controller cannot be 
readily tuned to each propeller, the comparison of thruster 
units will be made with the data from the torque mode tests. 
Ideally, the steady state power consumption of the velocity 
mode can be tuned to a value close to the power used in the 
torque mode. The deficiency in this controller will be 
commented on in Chapter 8. 

In general, the MFM motor with the 10:1 gearbox surpasses 
Hic Pittman motor «with, “the Arl- gear- reduction. The only 
exception is the JASON propeller, which performed best with 


the Pittman motor unit. Since the motors are electrically 
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comparable, it can be inferred that the 10:1 gearbox is a 
better choice for the 88, 86 and 316 propellers. 

The JASON propeller has the most limited range of thrust, 
regardless of power input. The maximum thrust for CHi 
propeller is near +/- 10 lbf. This gives insufficient thrust 
for ABE's sprint capability of 2 knots (requiring about 13 
lbf). Due to this limitation, the JASON Propeller will not be 
considered further. 

Figure 6.5 shows the comparison of the propellers using 
the MFM motor in torque control and the 10:1 gearbox. For the 
forward thrust direction below 13 lbf, the 88, 86, and sae 
propellers have nearly the same thrust/power characteristic. 
Since this is the region of interest, there is little to 
differentiate between the propellers. However, ABE will be 
Operating in an area of current gradients. This means that 
forward and reverse thrusts will be necessary to maintain 
constant velocity while traveling on a closed circuit 
trajectory. Reverse thrust will also be necessary while ABE 
1s maneuvering at docking. Therefore, the main propulsion 
thrusters should have the best possible reverse thrust 
efficiency without severely effecting the forward thrust 
characteristic. The 316 propeller has the lowest power 
consumption for any astern thrust demand. For comparison 
purposes, Figure 6.6 shows the velocity mode control curves 


for the MFM motor and all propellers. The tuning of the 
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controller increases the separation of the 316 curve from the 


other curves. 
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Figure 6.5 Propeller Comparison, Torque 
Mode, MFM Motor 


Velocity Mode Propeller Curves 


Input Power, Watts 


Thrust, lbf 





Figure 6.6 Propeller Comparison, Velocity 
Mođe, MFM Motor 
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6.4 Summary 


For steady state operation, the experimental thruster 
consumes less power than any other unit considered. In the 
range of - 10 lbf to + 10 lbf thrust demand, the experimental 
thruster has a power consumption advantage of 12% over the 
nearest competitor. This assumes that the thruster operates 
over a given trajectory, where any changes in thrust occur 
slowly and there is an equal demand of forward and astern 
thrusts. Chapter 7 examines the case where these assumptions 


are not valid. 
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Chapter 7 


Balancing Dynamics and Efficiency 


7.1 Introduction 


Chapter 6 compared the combinations of propellers and 
motor/gearbox units in steady state. Although this gives 
important insight into the performance of the thrusters, this 
1s only a small portion of the story. ABE will rarely operate 
ina "steady state" environment. During ABE’s flight along a 
preprogrammed trajectory, there will be accelerations, 
constant velocity runs and unknown current gradients. These 
factors will drive ABE away from the steady state toward a 
richly dynamic environment. 

In this chapter, the experimental propeller will be 
compared with the two model airplane propellers presented in 
Chapter 6. The experimental propeller will be driven by the 
MFM motor with the 10:1 gearbox. The airplane propellers will 
be combined with each of the two motor/gearbox units. The 
comparison will be based on the model of Chapter 4, with the 
parameters experimentally determined for each thruster. These 
parameters are shown in Table 7.1 . The following definitions 


apply Eo “Table 7.1 
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Thruster 1 is the experimental thruster consisting of the 
MFM motor, the 10:1 gearbox, and the 316 propeller. 
Thruster 2 is the MFM motor, the 10:1 gearbox, and the 88 
propeller. 

Thruster 3 is the same as thruster 2 except with the 86 
propeller. 

Thruster 4 and 5 are the same as 2 and 3 (respectively) 


except with the Pittman motor and the 4:1 gearbox. 


PU O i 


1 (M316) 


2 (M88) 


4 (P88) 


5 (P86) 





The response of each thruster was simulated on each of 


six trajectories: Hover, Track #1 through Track #5. The 
tracks (#1 through #5) are shown in Figure 7.1. In the Hover 
mode, the simulated vehicle is initially displaced 1 
centimeter off the desired position. The vehicle then 


attempts to regain and maintain the desired position. 
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Figure 7.1 Tracks #1-5, 


current. 
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The qualitative description of the tracks is as 


Track 
Track 


Track 


Track 


Track 


T: 


2: 


Low Velocity, 
High Velocity, 


High Velocity, 


High Velocity, 


Low Velocity, 


Low Acceleration, 
Low Acceleration, 


Low Acceleration, 


High Acceleration, 


High Acceleration, 
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Track <s:*hnas: vat. Sms oppostug 


EOL Lows: 


No Current. 
NO Current. 
Opposing 
Current. 

NO Currenecs 


No Current. 


7.2 Dynamic Comparison 


From step response data for each thruster unit under 
consideration, the speed of response Can be rated. The 
criteria is the time until the thrust output stays within 5% 
of the final value. The thrusters were judged for a thrust 


step of 50 Newtons (approximately 11 lbf). The results are 


shown in Table 7.2 


M316° 0785 "sec 





Table 7.2 Time Response " Experimental Thruster 


The 18-inch airplane propeller with an 8-inch pitch has 
the fastest response of all the propellers. For small thrust 
steps the differences between the speeds of response is large. 
For larger steps the time responses approach the same limit. 
The experimental propeller has the slowest response of the 


tested thrusters. 
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7.3 An Open Loop Force Controller 

In order to evaluate the thrusters a simple force 
controller will be implemented. This will be a PD controller 
with poles placed such that the system’s natural frequency is 


@, = 1.0 rad/sec, and the damping ratio is N = 0.707 
7.3.1 The Simulated Vehicle 


In order to test the thrusters on the aforementioned 
tracks, some vehicle dynamics must be considered. The vehicle 
modelled will be ABE. The drag of the vehicle is the drag 
calculated in Chapter 2. The drag at 1 knot 1s used to 


calculate an effective drag coefficient C', such that: 


Kg 


/ 1 Bad 
= PAC, = a (7.1) 


Cp 


The effective mass of the vehicle will be the 
displacement of the three main body sections of ABE plus an 
added mass term. The added mass is used to account for the 
volume of water in and around the vehicle that must be 
accelerated in order for ABE to move through the water. For 
this analysis, the added mass is assumed to come from a sphere 
of water in front of each body section. Each sphere has a 
diameter equal to the maximum diameter of the associated 


section. The displacement of ABE is approximately 800 kg and 
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the added mass is about 160 kg. This leads to an effective 
mass of 960 kg. 

Since three thrusters are used to drive the vehicle, and 
the point on their operating curves is significant to the 
analysis, it is assumed that the three thrusters share the 
load equally at all times. This allows the simulation to be 
run for one thruster using one third of the values calculated 
above. Therefore each thruster will ‘see’ an effective drag 
coefficient of C’, = 22.67 kg/m and an effective mass of 320 


ge 


7.3.2 Power Consumption Comparison 


Each thruster was run on each simulated track. The 
average power during the run was computed and tabulated in 
Table 7.3. Initially, it 1s assumed that the thruster has the 
same forward and reverse power to thrust characteristic. 
Double line blocks indicate the lowest mean power consumption 


for a particular trajectory. 


90 


EEES 


Track 
| 1 (M316) 19.96 | 24.56 DO ila fuel 86 [oss | 59 


2 (M88) EE al 49596 1 156 ii y 27:24 | 
3 (M86) 20619 Ee Ae PS 39.99 SS 28.445 


4 (P88) Pi 85 LIO DEI 47.98 40.50 LS 


5 (P86) 2 O 25.00 A 40.99 SA 25-96 


Table 7.3 Mean Power Consumption in Watts 

















The assumption of equal ahead and astern performance is 
not valid (see Chapter 6, Figure 6.5). A penalty factor must 
be applied to each thruster to correct for the increased power 
consumption in the astern thrust direction. The penalty 
factor (PF) used in this analysis is the ratio of slopes from 
the steady state power versus thrust characteristic for each 
thruster (Figure 6.5). The penalty factor for each thruster is 
shown in Table 7.4 . In order to apply this correction, some 
assumption must be made about the ratio of time at ahead 
thrust versus the time at astern thrust. When it 1s assumed 
that the thruster operates between 1/4 and 1/2 of the time 
astern, the results for Tracks #1,2,4,and 5 are the same in 


PCrMS OL meNrUStCers-rankwings (1er “Ene “sensitivity Eo. this 
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assumption is low) Except as specifically noted, it is 
assumed for a given track, that one third of the time is spent 
under astern thrust conditions. It is further assumed that 
the astern thrust distribution is the same as the ahead thrust 
distribution. Track 3 and Hover mode are the exceptions. 
Since Track 3 has a constant opposing current, it is assumed 
that the thruster only operates in the ahead direction. For 
the Hover mode, it is assumed that the thruster operates for 
equal times at ahead ana astern thrust. 

In addition to the correction for the asymmetric nature 
of the propellers, the quiescent power load (PQ) of the 
controller must be removed. This load is the power the 
controller draws in spite of the load from the motor 
operation. For the controller used in these tests, this load 
list PO: 12922 watts. This 19.2 Watts is high ana obscures 
the effect of the propellers. The equation used to apply both 


torrecilons is: 


A (VITO) a (70920 
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t (M316) 










2 (M88) 


sd 


1.4 Penalty Factors 





In order to evaluate the power used by the propeller and 
controller to accelerate the vehicle mass, it is helpful to 
remove the portion of the mean power used to overcome the 
vehicle drag. This mean power is calculated by integrating 
(Drag Force)* (Velocity) for each track, and then dividing by 
the simulation duration. The resulting mean power data after 
correcting for asymmetry, and removing the quiescent and 
'drag’ power, 15 shown in Table 7.5. The last row of this 
table is the calculated mean power required to overcome the 
drag for each track. Double boxes indicate the lowest mean 


power for a given track. 
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Tr 11.91 Ter 17.77 || 10.66 
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AAA 
2.14 9.43 SITU: 













Table 7.5 Corrected Mean Power Data in Watts 

Table 7.3 shows that the 88 propeller is tied for first 
place with the 316 propeller. When the 88 propeller is 
corrected, its loss of efficiency in reverse thrust operation 
drops this propeller to second place. This conclusion holds 
for several correction formula weights between 1/4 and 1/2. 
All the weights calculated gave results similar to those 
presented in Table 7.5. 

Figure 7.2 shows a break down of mean power consumption 
for each track and Ehalsvers i 2 fondeoe The cross-hatched 
section up to 19.2 Watts indicate the controller’s quiescent 
power load. The solid filled region immediately above the 


19.2 level indicates the power required on each track to 
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Mean Power Consumption 


Thruster Unit 





Figure 7.2 Power Consumption 


overcome vehicle drag. The enclosed area remaining above this 
level accounts for the remaining power used by the controller 
and the propeller to accelerate the vehicle and overcome the 


hydrodynamic losses. 


7.3.3 Comparing the Dynamics 


The 316 propeller has the best power consumption over the 


tracks tested, with the exception of Hover. The Track 2 


velocity response of the 316 propeller is shown in Figure 7.3. 


Figure 7.4 shows the same response for the 88 propeller. 
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The position response for both propellers in Hover mode is 


shown in Figure OF 


M316 Propeller, Track #2 


Time, Seconds 


Figure 7.3 Experimental Thruster Track 
#2, Velocity Response 





M88 Propeller, Track #2 


15 20 25 





Time, Seconds 


Figure 7.4 The M88 Thruster, Track #2, 
Velocity Response 
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Hover Position Response 
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Figure 7.5 Thruster Hover Response 





7.4 Summary 

The simulations presented in this chapter show there is 
a trade-off between power conservation and rapid dynamic 
response. This is clearly illustrated by the data for the 
Track 2 simulations presented in Table 7.5, Figures 7.3 and 
7.4. The M88 thruster follows the velocity track better than 
the M316 thruster. The «cost OL this “slightly improved 
response is 2 Watts, a 14% increase over the mean power 
consumed by the M316 thruster. This 2 watts, when multiplied 
by three forward propulsion thrusters, is 6% of the total 
propulsion budget for ABE. 

Better dynamics is a more important criteria when 


evaluating the hover mode thrusters. By responding quickly, 
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the M88 thruster reduces the magnitude of the hover limit 
cycle. This in turn limits the vehicle velocity and power is 
saved. This is shown by the plot in Figure 7.5, where the M88 
thruster uses 40% of the power used by the M316 thruster. 
For a typical ABE track having long constant velocity 
legs, rapid dynamics are less important. For a well developed 
controller, providing slow acceleration up to survey speed, 
the M316 thruster is the optimum choice. The cost in dynamic 
response is negligible in this case and the power savings is 


significant (14%. for Tracks rana 
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Chapter 8 


Summary, Conclusions, and Recommendations 


The purpose of this thesis was to develop an optimized 
propulsion system for ABE. From the beginning, special 
attention was given to maximizing the power efficiency of the 
system’s components. In this light, a propeller was designed 
specifically for the ABE vehicle and its predicted mission 
profile. This custom propeller was then matched to a suitable 
DC brushless motor. These units were connected via a speed 
reducing gearbox, thus taking advantage of the high efficiency 
operating regions of both components. 

Following its design, the ABE thruster was modelled using 
non-linear lumped parameter modelling techniques. The final 
model was reduced to three parameters (C,, o, B) describing 
the thruster. Using this model, five different thrusters were 
evaluated, including the designed ABE thruster. The three 
parameters for each thruster were determined experimentally, 
using system identification methods. From the experimental 
data and the modelled response, the static and dynamic 
performance of the thrusters was compared. Finally, a simple 


controller was implemented, and a simulated ABE vehicle with 
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model thrusters, was run on several tracks. This allowed 
evaluation of the control performance of each thruster. 

The result of this analysis indicates there is a trade- 
off between control performance and efficiency. In general, 
the designed propeller had the highest efficiency and the 
slowest response. The propeller with the fastest response, 
demanded more power in order to achieve this improved 
performance. The exception is in hover mode where the 
improved response limited the vehicle limit cycle and resulted 
in power savings. The simulated trajectories were typical of 
the predicted ABE track profiles. From these simulations, it 
became evident that the power used to overcome hydrodynamic 
drag is nearly insignificant for the proposed ABE operating 
velocities. The overwhelming share of the power goes first, 
to powering the control electronics, then to the un-modelled 
hydrodynamics and the power required to accelerate the 
vehicle's mass. Therefore, the use of the typical drag 
calculation is not sufficient to predict the propulsive power 
requirements for ABE. 

As the track length increases from the simulated track 
length, which is on the order of 10 meters, the importance of 
improved propeller dynamics is overshadowed by the need for 
higher efficiency. For ABE, this means that the designed high 
efficiency propeller is the best choice. The power cost for 
the improved propeller response is too high (6-14% as 


modelled) and the improvement in velocity tracking is small. 
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Additionally, the limitations imposed by the slower response 
of the designed thruster can be compensated for by the 
Pmeellagent “application. of control algorithms. Further 
research should start by determining the best way to control 
the thruster in order to best exploit its improved efficiency 
and to correct for the slow response. Following  clesely 
behind additional control research, should be the coordinated 
planning of the ABE trajectories. These tracks must be 
developed in a manner that allows the control system to 
operate the thrusters at peak efficiency. This includes using 
low accelerations and maximizing constant velocity legs as 
much as possible. Simulation runs with rapid acceleration to 
speed, consumed up to 60% more power than low accelerations to 


the same velocity. 
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